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Abstract 

A Mach-Zehnder interferometer with absorbing Fabry-Perot cavities is used to measure 
the optical phase differences upon rotation around a vertical axis. The details of the 
expected experimental results are described based on the idea that a conducting medium 
can introduce phase differences due to changes in material conductivity that are expected 
to occur when a medium moves with respect to a reference frame. Experimental details 
and results are discussed. 


1 Introduction 

In Einstein’s theory of special relativity |1] it is assumed that light propagates in empty 
space with a constant velocity independent of the velocity of the observer. Many experi¬ 
ments support this assumption. However, it is well know that most of these experiments 
are also supportive of the assumption that light propagates through a medium under 
conditions of time dilatation and Lorentz contraction. That light, or in general electro¬ 
magnetic waves, can be interpreted as waves propagating through a medium is the basis of 
the success of the Maxwell equations. Lorentz [2] was able to show that Maxwell equations 
could lead to Lorentz contraction and time dilatation for systems moving with respect 
to the medium. When the medium is referred to as ether , this is known as the Lorentz 
ether theory. Until recently is was thought that Einstein’s theory of special relativity 
and Lorentz ether theory were experimentally indistinguishable. The preferred reference 
frame or ether frame could never be detected and hence ignored. However, it has become 
clear that in principle experiments are possible to distinguish between them 0.0 

One type of such experiments is based on the occurrence of Wigner rotation or Thomas 
precession due to the fact that two successive non-co-linear Lorentz transformations can 
only be described by means of a Lorentz transformation of the composite velocity and a 
rotation over the Wigner rotation angle [5]. This in principle should enable the detection 
of the ether frame. 

Another type is based on experiments where light has some interaction with matter. 
As soon as a medium interacting with electro-magnetic waves exists, it is possible that this 
interaction enables the detection of the ether frame. A theoretical possibility has recently 
been outlined by Spavieri [6] and Consoli (7] although the first one has been disproved by 
experiment [8]. Another possibility is the occurrence of superluminal signal transport, as 
assumed to be possible in quantum mechanics due to its non-local character as discussed 
for instance by Einstein [9] for the Einstein-Podolsky-Rosen thought experiment. 

The absorption process of light inside a conductor is poorly understood. It is based 
on Ohm’s law relating the free current in a conductor with the electric held strength, 
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interpreted as a resistance against the vibrations of the electrons. In classical theory this 
resistance can be interpreted as due to scattering of electrons at certain time intervals, 
so that the electron momentum and energy is transferred to the molecules, ions or atoms 
constituting the medium. It can also be interpreted as a frictional force due to interaction 
of the electron with the atom to which it is bounded. In current theories quantum 
mechanics is involved to find models for the resistivity. However, none of these models are 
proven to be explicitly Lorentz co-variant and hence conductivity is a material property 
that could depend on the velocity of the material with respect to the ether. 

According to some authors [10], [IT], pj2|, |T3] the interaction between light and matter 
results in super luminal signal transport. Longhi ra shows that light transmitted through 
two fiber Bragg gratings at some distance apart experiences a phase shift independent of 
this distance. Olkhovsky [14] argues that this also holds for resonant tunneling through 
a double barrier. Each barrier could be build from a fiber Bragg grating as was done by 
Longhi. In such a case the wavelength must be carefully tweaked so that the transmission 
is small (but not too small). However, any double barrier system (for light or quantum- 
mechanics) can experience such a phase behavior as long as the total transmission is small. 
In principle there is no reason why the barriers should consist of fiber Bragg gratings. A 
simple metal layer on a glass substrate (appropriately chosen thickness) gives similar 
results as has been shown by for instance Monzon [15], Giust [16j and de Haan S3. 

The above arguments can be combined and resulted in the experiment described below. 


2 Refractive index of a conducting material 


In a simple model of a layer made of a conductor with a conductivity a, the plane wave 
solution of Maxwell’s equations can be developed from the Lorentz-Drude model [18], 
where the relation between the polarization of the material, P(f, t ) is proportional to the 
distance of each electron Xi to its position when the electric held E(f,t ) would be zero. 
This gives 

P(r, t ) = n e qx(r, t ) 

where n e is the electron density of the material, q the electric charge and x(r, t ) the 
average distance to the equilibrium position of all electrons. Here, it is assumed that 
the number of electrons in a region where P(r, t ) is (almost) constant is large enough 
to enable a statistical averaging with small enough standard deviation. The equation of 
motion of each electron is assumed to be (ignoring the magnetic held and the magnetic 
interactions) 


+ Kx(f, t) + f = , (B(r, t ) + -?{r, t) 


dt 2 


dt 


where e G is the permittivity of vacuum, m e is the electron mass and k e is an elastic constant 
of an elastic force driving the electron back to its equilibrium position, a is due to the 
interaction with all other electrons and ions in the material and for a spherical cavity in 
a homogeneous polarized material equal to approximately 1/3. The friction coefficient / 
can be interpreted as due to the collisions of the electrons with the atoms. If the mean 
time between collisions is r then / = m e /r. 

A conduction current J = qn e v can be interpreted as a how of free electrons (a = 
k e /(q 2 n e )) moving with constant uniform speed v so that the acceleration is zero, hence 
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the above reduces to fv = qE. According to Ohm’s law J = crE so that / = q 2 n e /a. 
Note that this interpretation does not hold for conditions where the electrons are not free, 
as then the friction coefficient might be small although the global conduction is negligible. 
In the harmonic case the above equation reduces to 

(-w 2 /(o;|e 0 ) - a!/e Q + iu /a) P(r,t) = E(r,u) 


where u p is called the plasma frequency equal to q\Jn e / ( m e e 0 ) and a' = a — k e /(q 2 n e ). 
Under consideration that D(r,u) = e 0 E(r,u) + P(f,u) and the constitutive relation 
between the electric held and the electric displacement D(r,u ) = e r (u)e 0 E(r,u) gives the 
relative permittivity of the material as 


e r (tu) — 1 + 


—u 2 /u 2 — a' + inu/u p 


where n = u p e 0 /a = l/(rcOp). 

The refractive index of the material is given by 


n(u) = \J p r e r (u ) = n r + in* 


where /i r is the relative permeability of the material, n r the real part of the refractive 
index and n, : its imaginary part. As before, the magnetic interactions are ignored so that 
pL r = 1. For a good conductor both a' « 1 and k < 1 so that the refractive index can 
be approximated by 
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for u r < 1 and by 


Hi = 



sju 2 - 1 


u r 


K 


2u 2 sju 2 - 1 


for u r > 1. Note that for u r < 1 the real part of the refractive index is proportional to k, 
hence inversely proportional to the conductivity. 

The real and imaginary refractive index for bulk Silver |19].[20] as function of wave¬ 
length (in vacuum) is shown in figure |T] For longer wavelength the above model of the 
refractive index is fitted to these data and shown as the solid lines. The plasma frequency 
u p = 1.2 x 10 16 1/s, or 7.8 eV corresponding to a wavelength of 0.15 /iin. This is similar 
to the value reported by Rakic [21 1 (9.0 eV), but he used a far more extensive model. 
Further a' = 3 x 10~ 5 and the fitted conductivity a = 1.3 x 10' S/m so that k = 8 x 10” 3 . 
The fitted conductivity is smaller than the DC value of 6.30 x 10 7 S/m at 20°C [22]. The 
bound electrons that are capable of oscillations experience a larger frictional force than 
the free electrons responsible for the DC conduction. For smaller wavelength the model 
is not correct any more. It could be improved by incorporating more oscillators (as has 
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Wavelength in micrometer 


Figure 1: Real (black) and imaginary (red) part of the refractive index of Silver as function 
of wavelength as reported by [T9] (circles) and [20j (squares) compared to the simple 
Lorentz-Drude model as described in the text. The inset is the same with a logarithmic 
scale for the refractive index to make clear that for short wavelength the model collapses. 


been done for instance by [21]). The wavelength dependence of the refractive index for 
some other metals with a good conductivity (such as Copper, Gold and Aluminum) is 
comparable to the one of Silver. For the visible light region it can be fitted to the same 
model, yielding a slightly different plasma frequency. 

3 Absorption in a conducting layer 

The general solution for a propagating plane wave in a homogeneous medium is 

T(r, t ) = cos [ojt ±k-r + (j)j 

where r is the position vector, t is time, T(r, t) represents any component of the electric 
or magnetic held, T, the amplitude of the plane wave at the origin, 0, the phase of the 
wave at the origin, oj the frequency of the wave and k the wave vector in the material. 
The ± indicates the direction of the plane wave, — in the direction of k and + against it. 
The dispersion relation is given by the real part of the refractive index 

kc 

— = n r 

LU 

where c equals the velocity of light in vacuum. The extinction coefficient a is given by 
the ratio between the imaginary and real parts 
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Figure 2: Wavelength dependent light-transmission measurements of several Silver lay¬ 
ers on top of a BK7 glass substrate (circles) compared to simulations for several layer 
thicknesses (colored lines). The black line represents the results for a 25 nm silver layer 
calculated with the simple dispersion model as described in the text. 



The electro-magnetic interaction of light with inhomogeneous materials can be calculated 
by division in regions where the material is homogeneous, solving Maxwell’s equations in 
each region and connect the solutions at the boundaries ram In this way, the reflec¬ 
tion, absorption and transmission of any optical system can be calculated and compared 
to experiment. 

For instance, wavelength dependent measurements of the transmission of a thin silver 
layer on top of a 3 mm thick BK7-glass substrate are shown in figure [2] and compared 
with calculations for several layer thickness using the index of refraction as given by the 
data points in figure [D It is clear that the transmission depends very strongly on the 
wavelength. Qualitatively, the transmission corresponds to that of a Silver layer with 
a thickness between 20 and 30 nm on top of a glass substrate of 3 mm. The black line 
represents the calculations for the simple dispersion model as described before and a Silver 
layer thickness of 25 nm. Again for smaller wavelength calculations differ considerably 
from the measurements. This is presumably because the used model refractive index 
differs from the actual one. Deviations in the refractive index can occur because the 
Silver layer is far from perfect. The layer tends to form small islands. A scanning electron 
microscope picture of one of the layers is shown in figure [3] The islands are clearly visible 
as the bright parts. Here, it is assumed that the effect of the islands can be incorporated 
by an effective thickness and refractive index averaged over the surface of the layers. As 
can be seen from figure [2] these parameters will change from sample to sample. In the 
following it is assumed that the refractive index of Silver is given by the Lorentz-Drude 
model with the fitted parameters mentioned above. 

Let 0 be the phase a light beam acquires when traversing a length L through the 
material. Then, the phase difference between phases acquired by a wave traversing a 
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Figure 3: Scanning electron microscope picture of a 25 urn Silver layer on top of a BK7 
glass substrate. The silver layer is created by islands visible as the bright parts. 


length L through the material when the conductivity changes with A a is 

Act 


A (j> = S a 


a 


where 


d(f) 

Sa = cr— 
da 


is called the conductivity sensitivity of the light beam phase to a relative change in the 
conductivity. A similar plasma frequency sensitivity can be introduced to indicate the 
change of phase when the plasma frequency is changed 


Sr. t, rl UJr, 


dcj) 

du n 


For a length of 25 nm through Silver for a wavelength of 0.633 pin S a = —n r kL = 
—0.016 rad and S u = 0.031. The plasma frequency sensitivity is larger than the conduc¬ 
tivity sensitivity due to the different dependency of the real part of the refractive index 
on a and co p . 

When passing through the Silver the beam of light is reduced in intensity, as the 
transmission is T = exp(—2nikL ) in this situation T = 0.146. By increasing the thick¬ 
ness the sensitivity increases but the beam intensity decreases exponential. The optimal 
performance depends on the minimal transmission for which the phase can still be de¬ 
termined with sufficient accuracy. Let us assume that this transmission is given by T min , 
then kL max = ln(l/T min )/(2ni) and = /n(l/T min )/(2a). Hence, for an optimal 

performance the ratio between the real and imaginary part of the refractive index should 
be a large as possible and as such the extinction coefficient as small as possible. One 
should keep in mind that this rule only applies for materials for which the main part of 
the refraction is determined by the dissipative process of conduction. In such a case the 
extinction coefficient is given by 


a = 


2oy 

K 


(1 - U>r) 



2a! \ 


for ay < 1 and by 


K 

2 - 1 ) 


for oj r > 1. For a frequency close to the plasma frequency (ay ss 1) the extinction 
coefficient changes drastically from very small values for frequencies above the plasma 
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frequency (a = 1.3k) to a maximum value at ay = l/\/3 (a = 0.77 /k). In the transition 
region around the plasma frequency (with a width of uj p k) the extinction coefficient has 
a discontinuity from very small to very large values for increasing frequency. However, 
for frequencies above the plasma frequency the real part of the refractive index does 
not depend on the conduction any more and hence there is no phase change when the 
conduction changes. Only when the frequency is below the plasma frequency the real part 
of the refractive index of the material depends on the conductivity. 

The transmission and sensitivities for a more complex model can be determined by 
numerical calculations following the method outlined in nn, For the sensitivity first, the 
phase of the transmitted wave is calculated for the model at hand, yielding 0i. Then, 
the conductivity of the model is increased by a small amount A a (about 1 ppm) and the 
phase is calculated again, yielding 0 2 . Then, the sensitivity can be calculated according to 
S c r = (02 — 0i)cr/Acr, yielding in this case S a = 0.062 rad. For a light beam traversing a 25 
nm Silver layer in air the transmission increases a bit 0.154 due to the interference between 
forward and backward traveling beams in the Silver layer. The conductivity sensitivity, 
S a changes sign and increases to 0.031, while the plasma frequency sensitivity in creases 
to a larger value of 0.81. The larger increase of the plasma frequency sensitivity with 
respect to the conductivity sensitivity is that due to the interference between the forward 
and backward light beams the imaginary part of the refractive index becomes important 
for the phase of the light beams. The relative influence of the plasma frequency change on 
the imaginary and real part is if the same order, however as the imaginary part is much 
larger, the absolute change in phase will also be much larger. 

Such a thin layer must be deposited on a glass substrate, influencing the interference 
further. Then the sensitivity and the the transmission will depend on the glass thickness. 
For a Silver layer of 25 nm thickness on BK7 glass of a thickness varying between 3 mm and 
3.004 mm (corresponding to a wavelength range in glass) the transmission and sensitivity 
is shown in figure [4] Depending on the glass thickness the transmission varies between 
0.14 and 0.3, the conductivity sensitivity, S a varies between 0.024 and 0.032 rad and the 
plasma frequency sensitivity between 0.7 and 1.2. For a real glass substrate the thickness 
will not be constant and hence the transmission and sensitivity must be averaged over 
the applicable thickness range. 

The effect can be enhanced by putting additional layers of Silver in series in the light 
beam. When two of these layers are put in series they constitute an absorbing Fabry-Perot 
cavity. 

4 Absorbing Fabry-Perot cavity 

With two semi-transparent mirrors a Fabry-Perot cavity can be constructed. Here, the 
semi-transparent mirrors are made out of a thin single layer of metal (25 nm of Silver) on 
a piece of glass (BK7) as schematically shown in figure 0 The mirrors can be positioned 
in four different ways to form the cavity: (A) glass substrates facing each other, (B) 
and (C) glass substrate facing metal layer and (D) metal layers facing each other. The 
length of the cavity can be varied. Normally the transmission of the cavity is maximized 
and the finesse should be very high. Here, the finesse should not be too high, as some 
light must be transmitted for the measurement. Hence, the metal layers are sufficiently 
thick to absorb part of the light and sufficiently thin not to be opal. Here, a thickness 
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Figure 4: Left: Transmission of light passing through a Silver layer of 25 nm on a BK7 
glass substrate as function of the relative glass thickness (this is the fractional part of the 
thickness expressed in wavelength). Right: Sensitivity due to variation of conduction (red 
line divided by 10) or plasma frequency (blue line) of the same model. 
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Figure 5: Schematic of Fabry-Perot cavity with length, L. Each semi-transparent mirror 
consists of a thin metal layer on top of a glass substrate. The mirrors can be positioned 
in four different ways to form the cavity: (A) glass substrates facing each other, (B) and 
(C) glass substrate facing metal layer and (D) metal layers facing each other. 
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Figure 6: Transmission (left) and phase (right) of light passing through a Fabry-Perot 
cavity A to D as shown in figure [5] as function of the relative cavity length (this is the 
fractional about of the cavity length expressed in wavelength). The blue line (No) gives 
the same quantities when the absorption is ignored. 


of approximately 25 nm was chosen and a wavelength of 0.633 /im, corresponding to the 
wavelength of a Helium-Neon laser. The results for the transmission and phase of a light 
beam passing through a cavity is shown in figure [H] for the situations A to D as described 
above. For each cavity the transmission has a maximum value when the standing wave 
matches the cavity. In that case the transmission is much larger than in the case for 
a single layer (values shown in figure [2] for 0.633 /irn). The minimum value is of the 
order of 1.5 % corresponding to the transmission through two single layers. The phase 
dependence shows that the phase of the light exiting the cavity behaves a-typical. When 
the transmission is low, the phase decreases when the cavity length is increased. When 
the transmission is high, the phase increases when the cavity length is increased. Note 
that the phase is relative to that of a light beam traveling the same distance through 
air. When the cavity would be made of non-absorbing materials only, the phase would 
be almost constant as shown by the blue lines. Almost, as also the multiple reflections 
inside the cavity influence the phase behavior, although to a much lesser extend. If also 
the reflection would be ignored the transmission becomes 1 and the phase constant. 

When for a fixed cavity length the conductivity of the Silver layers changes, also the 
transmission and the phase change. The sensitivity, S a as function of the cavity length 
is shown in the left graph of figure [7] for cavity types A to D. The maximum sensitivity 
occurs when the transmission is lowest. When somehow the conduction is constant, but 
the plasma frequency changes, the sensitivity is shown in the right graph of figure [TJ 

5 Mach-Zehnder interferometer 

If such Fabry-Perot cavities are put in each arm of an optical (wavelength, A) Mach- 
Zehnder interferometer (see figure [HD, then the output signal is determined by the intensity 
and phase difference between the two optical beams. 

The function of the beam splitters can be described by using the transfer matrix 
approach mm- Let A and i 2 be the amplitudes of the light waves at the two inputs of 
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Figure 7: Sensitivity, S a as function of the cavity length of Fabry-Perot cavity A to D 
as shown in figure [5] as function of the relative cavity length. Left: sensitivity due to 
conductivity changes; Right: sensitivity due to plasma frequency changes. 



Figure 8: Mach-Zehnder geometry for double Fabry-Perot cavity. 


the beam splitter and cq and o 2 the same for the outputs (see figure [§]), then for an ideal 
beam splitter 

As) 

The 2x2 matrix is referred to as transfer matrix and describes the action of the beam 
splitter. The incident light wave present at the first beam splitter, A is provided by a He- 
Ne laser with a wavelength of 632.8 nm. When only one light beam is applied, the output 
of both exit waveguides can be inferred by taking i\ — 1 and i 2 = 0. Then according to 
the above equation 

( 02 ) = ( -i ) 

The wave energy is split over both exits. The phase difference between the waves of 
both exits is — 7t/2. The propagation of the wave through the arms of the interferometer 
passing the Fabry-Perot cavity until the second beam splitter can be described by a 
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simple multiplication of the amplitude at the beginning of the arm by the transmission 
and a phase factor describing the optical path length experienced by the light. Here, 
it is assumed that the transmission of the first arm is T\ (complex to take into account 
the phase difference) and t 2 for the second arm. At the second beam splitter the waves 
interfere and result in the waves at the exits denoted by o 3 and 04 according to equation 
the above equation where i\ and *2 are replaced by i 3 and i 4 and similar replacement for 
oi and o 2 

03 V 1 ( 1 -i\ f h 

04 j a /2 V 1 / V *4 

If the above equations are combined one gets 

03 \ 1 f 1 W H 0 W 1 
o 4 J 2 { -i 1 J { 0 t 2 J { -i 


The sum of the intensities at the exit of the last beam splitter can be written as 


04 1 2 + |o 3 | 2 


T 1 + T 2 

2 


where Tf = |ti | 2 is the transmission of the Fabry-Perot cavity in arm 1 and a similar 
relation for arm 2. The transmission of both cavities can vary with the length of the 
cavity according to figure [ 6 ] depending of the type A to D. When the glass thickness 
varies one has to take that into account the by an appropriate averaging. The visibility 
of this interferometer can be defined as the relative difference between the two outputs of 
the last beam splitter, hence 
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where <p\ = argTi is the phase acquired by the light beam while passing through arm 1 
and a similar relation for arm 2. Here also, the phases (j) 1 and <f> 2 vary with the same 
dependence as the transmissions. 

To check if the above described effects are observable in an optic interferometer, mea¬ 
surements were performed with the Mach-Zehnder interferometer corresponding to figure[ 8 ] 
with Fabry-Perot cavities as described in the previous section. The light of a stabilized 
He-Ne laser (type Coherent 200, linear polarized, 0.5 mW, maximum mode sweep 10 
MHz) was coupled into the i\ arm of the interferometer via a polarization dependent 
optical isolator (isolation at least 35 dB). The interferometers were put in a tempera¬ 
ture controlled environment where the temperature control was within 3 mK. One side 
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Figure 9: Sketch of beam splitter geometry. 
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of each Fabry-Perot cavity was mounted by means of a piezo stack enabling the indepen¬ 
dent change of the length of the cavities. The light intensity at outputs 03 and 04 of the 
Mach-Zehnder interferometer were measured by two amplified silicon detectors. The sum 
of these intensities relative to the output intensity of the laser as function of the voltage 
applied to the piezo stacks is given in the right graph of figure [TUI For comparison also 
the ideal theoretical values are shown in the left graph. The occurrence of minima and 
maxima is the same, but the sharpness of the features in the measurements is much less 
than for the ideal theoretical calculations. This is due to the uncertainties on the exact 
sample parameters. Especially the glass thickness variations play an important role as is 
shown in the middle graph of figure [TOl This graph was calculated by using a variation in 
the glass thickness of ±20 nm which is comparable to the surface flatness of the substrate 
used (65 nm). The visibilities can be calculated from the intensities at outputs o 3 and o 4 
according to the above equation. They are shown in the right graph of figure |TT] as func¬ 
tion of the voltage applied to the piezo stacks. For comparison also the ideal theoretical 
values are shown in the left graph. Here also the occurrence of minima and maxima is 
the same, but the sharpness of the features in the measurements is much less than for 
the ideal theoretical calculations. The influence of the same glass thickness variations as 
before result in the middle graph of figure HU For increasing V\ voltage it seems that the 
minima and maxima shift a bit to larger values of voltage V 2 . During the measurements 
voltage Vi was set to a fixed value and V 2 was changed from its minimum to maximum 
value. A small drift in the interferometer path length of one arm with respect to the 
other will result in this small drift (a complete measurement took about 80 min). The 
measurements show that the transmission of the Fabry-Perot cavities behave as expected 
and that it is possible to measure the phase changes due to variation of the properties of 
the cavities. This behavior is exploited in the experiment that is described in the next 
section. 
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Figure 10: Theoretical ideal (Left), theoretical non-ideal (Middle) and measured (Right) 
sum of intensities of outputs of a Mach-Zehnder interferometer as function of the voltages 
applied to piezo stacks. 



6 Sun and Earth velocity components 

Let us define the Sun and Earth velocity components as measured by an observer on 
Earth. The rotation speed of the Earth around its own axis will be ignored. 
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Figure 11: Theoretical ideal (Left), theoretical non-ideal (Middle) and measured (Right) 
phase difference in radian of a Mach-Zehnder interferometer as function of the voltages 
applied to piezo stacks. 
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Figure 12: Definition of celestial coordinate systems. Left: Celestial view; Right: Earth 
view. See the text for the explanation of the symbols. 
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Figure fl2l shows the definition of the different coordinate systems and their connec¬ 
tions. Two views are presented. The left view represents the celestial view with the Sun 
in the origin. The right view represents the Earth view with the Earth in the origin. 
The ecliptic or celestial coordinates are defined with respect to the Earth orbit around 
the Sun. A is the ecliptic or celestial longitude (red angle from vernal equinox). ft is the 
ecliptic of celestial latitude (red angle from plane through ecliptic). The so-called fixed, 
stars can be found in a direction for which (3 and A are constant. In vector notation the 
direction of the fixed stars is given in ecliptic coordinates by 



/ cos /3 cos A 
— cos /3 sin A 
\ sin ft 


The equatorial coordinates are defined with respect to the projection of the Earth equator. 
This projection makes an angle fitut with the ecliptic, cj> t n t is the Earth’s axial tilt with 
respect to the ecliptic. The rotation axis is the line between the vernal and autumnal 
equinox. In the coordinate system used, this corresponds to the x-axis. Currently 4>tut = 
23.44°, but the Earth axis processes slowly so that the equatorial coordinates only have 
meaning when the epoch of observation is given, a is the right ascension (blue angle 
along celestial equator) and S is the declination (blue angle going up from plane through 
celestial equator). In vector notation the direction of the fixed stars is given in ecliptic 
coordinates by 


/ 

Re = 

V 


cos 5 cos a 
— cos S sin a 
sin 5 


It is possible to transform the Equatorial coordinates into Celestial coordinates and vice 
versa by means of 
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With respect to an observer on Earth, the coordinate system used is shown in figure [T2J 
The Earth coordinate system is defined with respect to the local North. A is the Azimuth 
(black angle on Earth surface from North direction toward East, h is the elevation (black 
angle going up from Earth’s surface). In vector notation a direction in Earth coordinates 
is given by 


Ra 


( cos h cos A 
— cos h sin A 
\ sin h 


When the observer’s location on Earth is defined by its longitude (Ao) and latitude (fio), 
it is possible to transform the Earth coordinates into Equatorial coordinates and vice 
versa depending on the time (universal time tur) and date of observation by means of 


Re 


( — sin fio cos A sin A 
sin sin A cos A 
\ cos fio 0 
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Figure 13: Definition of the coordinate system of the observer at some position on Earth 
(Horizontal coordinate system). See the text for the explanation of the symbols. 
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where A = Aq + A Earth + 27t tur/Tcay is related to the sidereal time by 


tsid — T j 


Day- 


IA — 7T, 27t| 

~2tt 


where the function |a, b\ denotes the value of a modulo b. 

The ecliptic longitude of the Earth at observation time is given by 


A 


Earth 


= 2tt- 


AT< 


SE 


T, 


Orbit 


where To r ut — 365.26 days is the orbit time of the Earth around the Sun. A T$e is the 
time passed since the moment the Earth was at the latest vernal (or spring) equinox. It 
can be calculated by means of the Julian Date, JD by 


A T se = | JD + 343.85743, T< 


Orbit 


where JD is the count of days since noon Universal Time on January 1, 4713 BCE (on 
the Julian calendar). Hence, the Julian Date of January 1, 2000 12h UT equals 2451545. 
The Julian date can be inferred from the Julian calender (. DD/MM/YYYY ) by means 
of 


JD = DD + 30MM + 367YYYY +1721029 — 2Q + F 
where 


' 7(MM —2) \ | r (Q 


(3 P 

VT 


Q = YYYY + F 


(MM -3 


V 12 

and F(x) denotes the largest integer below x. 


and P = 1 + F 
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Parameter 

Symbol 

Unit 

CMB 

Miller 1 

Miller 2 

Ecliptical latitude 

(3Sun 

Degree 

171.9 

84 

-83 

Ecliptical longitude 

^Sun 

Degree 

-11.1 

159 

-41 

Right ascension 

C^Sun 

Degree 

168.2 

255 

73.5 

Declination 

3 Sun 

Degree 

-7.0 

68 

-70.6 

Magnitude velocity 

^Sun 

km/s 

369 

205 

205 


Table 1: Sun velocity parameters for several models. 


Let us define the velocity of the Sun through the Ether in ecliptical coordinates as 

( COS/3s un COS A Sun 
VSun,C VSun COS (3Sun sLl A Sun 

\ sin ft Sun 

where /3 s is the ecliptical latitude and A sun is the ecliptical longitude and vs un is the 
magnitude of Sun’s velocity through the ether. In the equatorial coordinate system this 
becomes 

l COS j3 S un COS A Sun 

VSun,E VSun sill 0tj.lt, Sm (3sun COS (fttilt COS (3sun sin A Sun 

V sin font cos f3 Sun sin \ Sun + cos <f.hut sin (3 Su n 
The velocity of the Earth with respect to the Sun in Ecliptical coordinated is given by 

/ Sin A Earth 


V Earth,C ^Earth 


COS A Earth 

0 


or in equatorial coordinates 


V Earth,E VEarth 


Sin XEarth 
COS (j) t ilt COS ^ Earth 
- sin dun cos ^ Earth 


where VEarth = 29.78 km/s. Lienee, the velocity of the observer with respect to the ether 
is v = vsun,c + vEarth,c, where the rotation of the Earth is ignored. This velocity can be 
transformed in the Earth coordinate system by applying the above transformation. 

The velocity of the Sun with respect to the ether is to be determined from the measure¬ 
ments. Miller [22] in his repetitions of the Michelson-Morley experiment found a preferred 
direction from his measurements, the values of which are shown in table [I] Miller 1 refers 
to his first derivation and Miller 2 to his more accurate determination in 1933. Both di¬ 
rections are nearly opposites. Another possible candidate is the reference frame in which 
the microwave background radiation is uniform (CMB) [26]. This direction is almost 
perpendicular to the first two mentioned. 


7 Effect 

Let us assume that the influence of the absolute velocity of the material with respect to 
the ether, v on the effective conductivity and plasma frequency is given by 
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where \r} a \ is of the order of 1. The value of q a depends on the way that the conductivity 
depends on the velocity with respect to the ether. For now this is unknown, but it should 
not differ too much from 1. Let (j) be the phase a light beam acquires when traversing 
through the Fabry-Perot cavities as defined in the previous sections. Then, the phase 
difference between phases acquired by a wave traversing when the direction of the light 
beam through the cavities changes from ( 0 to ( (both unit vectors) with respect to v is 


A (j) = S a r] a 


(C - Co) • V 
c 


When the interferometer is rotated along an axis represented by unit vector m, with an 
angle 9 the vector ( is given by the Rodriguez rotation formula 

C = Co + mx£sm0- (<fo - m(m • Co)) (1 - cos #) 


rotating the vector Co around a unit vector rh over an angle 6 according to the right hand 
rule. When („lm then this reduces to 

C = Co cos 9 + (± sin 9 


where (j_ = m x Co is in the plane of rotation perpendicular Co- Hence, 


A0 = (cos {9 — 9 a ) — cos 9 0 ) 

where j3 — \v x m\/c the magnitude of the component of the ether velocity vector in the 
plane of rotation relative to the speed of light and cos 9 0 = Co • v/\v x m\, the length of 
that component parallel to Co- Hence, for a complete rotation A0 changes according to 
a cosine with an amplitude of j3. The maximum occurs when 9 = 9 0 and the minimum 
in the opposite direction. Both (5 and 9 a depend on the magnitude and direction of the 
ether velocity with respect to the rotation plane and the initial direction of the setup. 

The normal to the rotation plane can be expressed in the Earth’s coordinate system 
according to 

/ 0 


m A = 


so that in the equatorial coordinate system 


V 



( 

— sin <t> 0 cos A 

sin A cos (j>o cos A \ 

( 

°\ 

/ cos (fro cos A \ 

rh E = 


sin (j) 0 sin A 

cos A — cos 4>o sin A 


° = 

— cos 4>o sin A 


V 

cos (j)o 

0 sin (j) 0 J 

V 

iy 

\ sin (f)o J 


Let us assume that Co is in the South-North direction so that its Azimuth is 0, hence in 
Earth coordinate system 

/ 1 


Co,A — 


so that in the equatorial coordinate system 


0 

VO 


Co,E ~ 


( — sin C>o cos A 
sin (j>o sin A 
\ cos 4> 0 
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Figure 14: Ether effect for an observer at a latitude of 51.8° and a longitude of 4.6° and a 
Sun velocity according to Miller 2 in table Q] as function of sidereal time (horizontal scale) 
and days after 1 January 2015 (vertical scale). Left: In-plane-of-interferometer ether 
velocity component (in 0.001 c); Right: Azimuth of maximum of effect upon rotation (in 
degree). 

From these equations the maximum magnitude and the direction of the maximum can be 
determined. 

For an observer at a latitude of 51.8° and a longitude of 4.6° (i.e. location Puttershoek 
in the Netherlands) and a Sun velocity according to Miller 2 in table [U the magnitude of 
the Ether velocity component in the plane of rotation of a horizontal interferometer (/3) 
and the Azimuth at which this maximum occurs are shown in figure HU It is clear that the 
effect changes depending on the sidereal time (horizontal scale) and season (vertical scale) 
as represented by the day after January 1, 2015. Both the magnitude of the maximum 
and the direction in which the maximum occurs change. 


8 Experiment 

The experimental set-up used is shown in figure [Q)J It is based on the assumption that the 
effect is first-order in the velocity with respect to the ether. When the light beams travel 
in opposite direction also the effect is opposite. When the Mach-Zehnder interferometer 
is constructed as shown in the figure then the effect doubles. The light of a laser (a 
stabilized He-Ne laser type Coherent 200, linear polarized, 0.5 mW, maximum mode sweep 
of 10 MHz) is coupled into a fiber optical cable (Patch Cable, FC/PC, 630 nm, PM Panda 
Style, 1 m, Thorlabs) by means of hberport F\ (FiberPort PAF-X-5-B, FC/PC & APC, 
EFL=4.6 mm, 600-1050 nm, Thorlabs). The fiber optical cable is led into a temperature 
stabilized containment schematical denoted by the dashed box. The temperature of this 
box is stable for days within 3 mK. For shorter times the temperature is stable within 
1 mK. The heat applied to control the temperature of the box varies between 1 and 4 W. 
This ensures that the possible temperature gradients within the inner compartment are 
less than 1 mK/m. Inside the box the laser light exits another hberport F 2 (similar to 
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50 cm 



F\). To prevent optical feedback from the interferometer into the laser a polarization 
dependent optical isolator O (IO-2D-633-VLP from Thorlabs, isolation at least 35 dB) 
was included. After transmission through the optical isolator the light beam is split into 
two arms at beam splitter B\ (Cube Mounted Polarization Insensitive Beam splitter, 400- 
700 nm, CM1BS013 from Thorlabs). The first arm (green line in the figure) is directed 
to a light detector Di (all detectors used are amplified Silicon detectors, PDA36A from 
Thorlabs). This detector is used to monitor the output intensity of the laser. The other 
arm continues to a second beam splitter B 2 where the light beam is split in arms A (red 
line) and B (blue line) of the Mach-Zehnder interferometer. These light beams are joined 
together at the third beam splitter B 3 after which the light beam interference is detected 
by detectors D 2 and D 3 as described in section [5l The light in arm A passes a Fabry- 
Perot cavity Cj similar to the one described in section |4l reflects at mirror Mi (12.7 mm 
diameter Protected Silver Mirror, 6.0 mm thick, Thorlabs PF0503P01) which is connected 
to a piezo crystal (Piezoelectric Actuator, Max Displacement 9.1 /irn, 3.5 x 4.5 x 10 mm, 
AE0203D08F Thorlabs) to be able to adjust the phase difference of the light at beam 
splitter B 3 . The light in arm B reflects at mirror M 2 (Cube Mounted Protected Silver 
Turning Mirror, Thorlabs CM1P01) passes a Fabry-Perot cavity C 2 similar as before, 
reflects at mirror M 3 (similar to Mi) into the direction of beam splitter B 3 where it 
is joined to the light of arm A. One window of each Fabry-Perot cavity is mounted by 
means of a piezo crystal (as before). All mirrors and windows are mounted with kinematic 
mounts to enable angle adjustments to make sure that the beams interfere in the cavities 
and at beam splitter B 3 . The set-up is designed in such a way that the lengths of arm 
A and B are approximately equal (both 464 mm within a few mm). This limits the 
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influence of laser instabilities on the measured phase difference. The phase is measured 
with an accuracy of approximately 0.03 rad. The influence of the temperature of the box 
on the phase was measured by means of an oscillating temperature (with an amplitude 
of 0.04 K and a period of 2.7 hours) and was 65(1) rad/K. This explains the need for 
a stable temperature. The transmission of both cavities is kept minimal by adjustments 
of the cavity length by means of the piezo crystals. In this way the transmission is kept 
minimal within 0.3 % of the transmission. 

For determination of the effect, the set-up is rotated along a vertical axis. Every 
hour the set-up rotates from 0 to -180 degrees, from -180 to +180 and from +180 to 
0 with steps of 30 degrees. For 0 degrees the interferometer is directed according to 
figure M For 90 degrees the arrow in this figure points to the local West. At every 
step the phase difference and transmission of the interferometer are measured during 20 
seconds to average out statistical variations. A complete rotation takes about 20 minutes. 
A typical response of the interferometer as function of angle is shown in figure Q+J The 
top graph of this figure shows the angle of the set-up as function of time. The middle 
one shows the phase difference of the interferometer during the same period as the upper 
graph. The bottom graph shows the transmission of the interferometer again during the 
same period. In the lower graph the phase oscillations due to the rotation of the set-up 
are clearly visible. 


9 Results 

At first, it seems that this interferometer gives an indication that there might be some 
effect on the phase difference due to the motion of the Earth, but the transmission is 
completely independent of the angle of the set-up. This corresponds to the anticipated 
effect. However, if such an effect exists it should depend on the time of day and year. 
Upon rotation of the set up, the amplitude and azimuth of the maximum should vary 
between certain minima and maxima depending on the orientation of the Earth velocity 
with respect to the preferred frame as discussed in the previous section. The Fourier 
transform of the phase difference (corrected for the linear assumed drift) as function of 
rotation angle gives this amplitude and azimuth of the maximum phase difference for all 
orders. The zeroth order is just the average phase difference during a rotation. The first 
order represents the amplitude and azimuth of that part of the signal that varies with 
the cosine of the angle of the set-up, corresponding to first order effects in v/c. The 
second order represents the amplitude and azimuth of that part of the signal that varies 
with the cosine of twice the angle of the set-up, corresponding to second order effects, 
and so on. The error in the values can be estimated from the difference between the 
Fourier transform of the data points measured for increasing set-up angles and the one 
measured for decreasing set-up angles (to find the systematic error due to the unknown 
drift) combined with the Fourier transform of the variances (to find an estimate of the 
statistical error). The amplitude and azimuth should vary with the sidereal time and 
epoch as discussed in the previous section. 

The results for 6 consecutive quarter years (Q1 to Q6) are shown in figure [T71 Al¬ 
though some variation is observed, the quality of the data is not sufficient to reach a final 
conclusion. The variations are less then are expected. This can have several causes as the 
absence of the effect or a masking effect or an wrong expectation of the magnitude of the 
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Time / minutes 


Figure 16: Top: angle of set-up as function of time. Middle: phase difference of in¬ 
terferometer during the same period as the upper graph. Bottom: Transmission of the 
interferometer again during the same period. 
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effect. 

It is possible to increase the sensitivity of the measurements by increasing the stability 
of the set-up (for instance mechanical stability and laser stability), reduce the rotation 
time of the set-up to reduce environmental effects and so on. It is recommended that 
the measurements are repeated at a higher altitude (as Miller has done) to reduce the 
influence of the possible entrainment of the ether that can reduce the effect considerably 
(note that the height above see level of Puttershoek is -10 m). 


10 Conclusion 

A novel experiment has been performed to observe possible deviations from the predictions 
of Einstein’s special relativity theory. Although the experiment reveals some sidereal 
deviations, the magnitude of the measured deviations is too small to reach any final 
conclusion. It is recommended that the experiment is repeated with a more stable set-up 
at higher altitudes. 
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Figure 17: Amplitude of phase difference upon rotation (left) and angle of maximum phase 
difference upon rotation (right) as function of sidereal time for 6 consecutive quarter years 
from April 2013 until November 2014. 
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